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Highly loaded turbomachinery tends to exhibit noticeable subsynchronous whirl 
at a frequency close to one of the rotor critical frequencies when the rotating 
speed is well in excess of that critical* The subsynchronous whirl is not a 
"forced** vibration phenomenon where the vibration frequency is the same as the 
forcing frequency; e.g. , rotor unbalancing. It is better classified as a posi- 
tive feedback phenomenon tending to amplify small disturbances at the system 
natural frequencies . The cause of whirl can be traced back to a number of forc- 
ing mechanisms. One of these is the variation of the blade tip forces due to 
nonuniform blade tip clearance. Alford (1) was the first to model this aerody- 
namic force on the blade. By assuming uniform upstream and downstream flow 
fields, he was able to argue that the variation of the blade tip force would be 
proportional to the efficiency variation. Therefore, if a turbine rotor has some 
small eccentricity (Fig. 1), then due to the efficiency variation, the blade 
force on the narrow-gap side must be greater than the blade force on the opposite 
side. The resultant force Fy on the rotor, acting perpendicular to the eccentri- 
city is then no zero. "This force can cause the rotor to whirl forward, i.e., in 
the same direction as the rotation. Qualitatively: 
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where 6t is the local tip clearance at azimuth e, e Q is the eccentricity and H 
is the blade height, D m is the mean diameter of the rotor and T n is torque on 
the rotor, 3 is a numerical factor of order unity depending on the blading. 


Contrary to Alford's assumptions, the presence of the blade tip clearance 
variation is bound to influence the whole flow field. The purpose of this paper 
is to investigate the effect of this non-uniform flow field on the Alford force 
calculation. The ideas used here are based on those developed by Horlock and 
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Greitzer (2). It will be shown that the non-uniformity of the flow field does 
contribute to the Alford force calculation. An attempt is also made to include 
the effect of whirl speed. The values predicted by the model are comnpared with 
those obtained experimentally by Urlicks(3) and Wohlrab (4). The possibility of 
using existing turbine tip loss correlations to predict £ is also exploited. 


The Model 

The method of approach is to first find the flow field due to the tip clearance 
variation, and then, using the velocity triangles, to calculate the resultant 
Alford force. The detail of the analysis is given in Ref. 5. It is to be stressed 
that a complete . three-dimens ioal investigation of the flow field with the 
non-uniform tip clearance is not attempted here. 

Consider a single-stage turbine with its rotor whirling at a constant eccen- 
tricity e Q (Fig. 2). If the hub/tip ratio is high, the flow field can be 
approximated to be two dimensional. Fig. 3 shows an actuator disc representing a 
two dimensional unwrapped blade row of circumference 2 ttR, where w is the relative 
flow velocity and c is the absolute flow velocity. If the Mach numbers upstream 
and downstream are small, the flow can be considered incompressible, as well as 
inviscid, outside of the blade row. 

The variation of the tip clearance around the annulus can be represented by 
6 | = E E n [i | n - Qti] (1) 


where Q is the whirling speed, y is the distance in the tangential direction, R 
is the mean blade radius, H the blade height, and the E n are constants. 


Since only the first harmonic is of interest, as will be explained later, 
higher harmonics are not included in the analysis and Ei = e Q /H. 

The perturbation velocities c x and c y are related to the stream 
function by 


C x 3y 


(2a) 



Upstream of the stag^, the flow is irrotational, 7 tyi = 0 and downstream the 
flow is rotational V $3 = The flow is steady in a reference frame 
rotating at the whirl speed Q, therefore the stream functions can be written as 

^ - A exp [ ^ (x + i y) - Q ti] upstream of the stage (3) 

1 K 
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= B exp [i ( - x + i y) - n ti] 

+ C exp [ i ( - x tan a i + y i - Q ti) ] downstream of the stage (4) 
R 


where 5.3 is the mean flow angle downstream of the rotor. Notice the 
irrotational term dies out far upstream and downstream. The vorticity downstream 
is convected by the mean flow at an angle 53. Again only first harmonics 
are considered.’ 

There are three matching conditions across the actuator disc surface: 

1) The axial mass flux is constant 


(5) 
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2) The relative flow leaving angle is constant 
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3) The third boundary condition involves matching perturbation quantities by 
using the known experimental information. 


The total- to-static turbine efficiency can be written as 
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where Cy 2 and Cy 3 are the y-components of the absolute flow velocity 
before and after the rotor (defined positive in the y direction), u is the blade 
speed, p 3 is the static pressure downstream of the disc, and poi is the 
total pressure upstream. 


Due to the local blade tip clearance variation, the efficiency can be written 
as the sum of the mean efficiency and its perturbation: 


n = rj + n 


(8) 


where to a first approximation 
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It is known that the second term in the above expression is small compared 
with the first term, if g is .not too far from the designed value. Hence the 
second term in Eq.(9) will be ignored in the analysis. 

We now perturb eq.(7) and use the additional condition that the stator leaving 
angle is also constant to eliminate c the third boundary condition is 
obtained. We also use the downstream y-momentum equation to eliminate the 
downstream pressure perturbation. In terms of the perturbation streamf unctions , 
the. three matching conditions can then be written as 
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Substituting the forms given by Eqs. (3) and (4) for i|>i and 3 , these 
conditions become three linear algebraic equations for the coefficients A, B and 
C. To obtain the blade force variation, the momentum equation is applied across 
the rotor blades. The tangential component of the force is given by 


f = p c ,(c - c ) H (11) 

y x3 y3 y2 

Since only the perturbation of force contributes to the resultant destabilizing 
force 


f =pc. (c.-c_)H+p, c 
y xl y3 y2 lx 

If we consider a reference frame 
together the forces perpendicular to 
net cross-force 


F = Real(f' e Qt ) cos (y/R) 
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( 12 ) 


rotating at the whirl speed 0, and add 
the instantaneous eccentricity, we find the 


d y 


(13) 
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This expression justifies the inclusion of only the first harmonic in the 
stream function expansion, since the higher harmonics contribute nothing to the 
Alford force. A non-dimensional excitation coefficient k X y is defined as 




eo 

H 


(14) 


where U s is the ideal total force in the circumferential direction on the 
blade, and e Q /H is the dimensionless eccentricity. 

Results 


Cases with stationary eccentricity are considered first. A typical plot of 
the axial velocity perturbation and the blade force variation vs. angular posi- 
tion are shown in Fig. 4. As was to be expected, these quantities are almost in 
counterphase with the tip clearance variation, which indicates that the nonuni- 
form flow field will tend to increase the Alford forces. 

Although Alford published his paper in 1965, very few experiments have been 
done to verify the proposed formula. Hence it is difficult to make any good 
comparison between the Alford model, the present analysis, and experimental 
results. . Two tests were done by Urlich (3) and Wolhrab (4), and more recently 
Vance and Laudadio (6) did verify the linearity between the Alford forces and 
eccentricity for a fan. Accurate experimental data seems very difficult to 
obtain due to the small magnitude of the forces to be measured. 

For the present analysis the data of Urlich seems to be the more suitable, 
since he also determined the efficiency vs. average clearance for an unshrouded 
turbine stage. From this information the values of A * 3r|/3(t/H) are found by a 
second order polynomial curve fit to the efficiency vs. clearance data (Table 1). 
Values of k^ are found by using A, together with appropriate flow angles (Ref. 

3) in the analytical model. The predicted values and experimentally obtained 
values of k^ are plotted in Figs. 5 and 6 , together with the values of k X y 
obtained from the Alford formula. It can be seen that the analysis gives a value 
of k^ higher than the Alford formula, indicating that the non-uniformity does 
increase the destabilizing forces. 

In addition, it is also clear from Figs. 5 and 6 that the new theory over- 
predicts the side forces, when compared to the experimental data. However, we 
should at this point also consider the effect of the nonuniform pattern of pres- 
sure acting on the turbine wheel rim. Fig. 7 shows that, although the pressure 
field just upstream of the rotor is in phase with the clearance variation, it is 
shifted by 90° just downstream, and that those nonuniform pressures act in such a 
way as to reduce the forward cross-force, 

A precise accounting of the effect of these pressure forces is not possible 
in the context of our actuator-disc theory. However, we can see that the 90° 
rotation of the pressure pattern must occur in the rotor passages, and so it is 
logical to calculate the induced side force by letting the rotated pressure wave 
act on an axial rim width equal to half a bladde chord. New excitation coeffici- 
ents, including this pressure force, are also shown in Fig. 6. Although the 
precise axial length over which the pressure forces act is difficult to deter- 
mine, it can be seen that they do contribute significantly to the net force, and 
that their effect is to essentially remove the overprediction. 
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The Use of Tip Clearance Correlations 


In many situations, curves of efficiency vs. clearance are not available. 
The possibility of using some of the existing tip clearance loss correlations is 
investigated here. The formula used is given by Dunham and Came (7), in the form 
of a tip loss coefficient which depends on flow angles and other stage variables, 
as well as on the clearance to chord ratio . This allows determination of n, and 
hence A for use in our model. Fig. 8 shows the predicted Alford forces based on 
the tip loss clearance correlation (Ref . (7) ) , together with the experimental 
results. It can be seen that the prediction is very reasonable, as it compares 
well with the data for small axial (stator-rotor) gap, although not as well for 
large axial gap. The fact that k^ varies with the gap between the rotor and 
stator is not well understood, and it may be highly dependent on the detailed 
geometrical arrangement of stator and rotor. 

Wohlrab (Ref. 4) tested a shrouded turbine. Since he did not supply data on 
efficiency vs. tip clearance, A has to be estimated, in order to compare the 
present analysis with the experimental values obtained. Encouraged by the re- 
sults for the unshrouded turbine, an estimated tip loss coefficient is used 
again, but this time the effective tip clearance is given by (Ref. (7)). 

-0.42 

t = (geometric tip clearance) x (number of seals) 

Since there is also a shroud, we would expect some seal forces acting on the 
rotor shroud band. These are found by applying a model developed by Lee, et al. 
(8). The results (Fig. 9) again show reasonable agreement, falling between the 
data obtained by Wohlrab for the two axial gaps tested. 

Effect of Whirl 


All the cases examined so far have been for an assumed static offset. The 
effect of whirl for a typical case is shown in Fig. 10. For a whirl speed equal 
to 50% of the rotating speed, the effect on k is small (± 10%). This is consis- 
tent with Wohlrab* s experiment, in which he found that within experimental accu- 
racy there is no difference between the values of k xy obtained from kinetic and 
static tests. It is noted, however, that the predicted effect of whirl (which 
could be expressed as an Alford damping coefficient), could not have been obtain- 
ed from the original Alford model, since it depends crucially on the delays 
between wheel motion and disc passage of the induced flow non-uniformities, and 
these were ignored in the Alford model. 

Conclusion 


The non-uniform flow field induced by the tip clearance variation tends to 
increase the resultant destabilizing force over and above what would be predicted 
on the basis of the local variation of efficiency. On the one hand, the pressure 
force due to the non-uniform inlet and exit pressure also plays a part even for 
unshrouded blades, and this counteracts the flow field effects, so that the 
simple Alford prediction remains a reasonable approximation. Once the efficiency 
variation with clearance is known, the present model gives a slightly overpre- 
dicted, but reasonably accurate destabilizing force. In the absence of effici- 
ency vs. clearance data, an empirical tip loss coefficient can be used to give a 
reasonable prediction of destabilizing force. To a first approximation, the 
whirl does have a damping effect, but only of small magnitude, and thus it can be 
ignored for some purposes. To gain more insight and understanding, more accurate 
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experimental determinations of the tip force must be made. It must be pointed 
out that the destabilizing force is also highly dependent on the geometry of the 
rotor and stator, as shown by the experimental dependence of the magnitude of 
destabilizing force on the axial gap between the stator and rotor. This effect 
of axial gap on the Alford force points to the need for more detailed study of 
the flow field around tip regions. 


References 


1. Alford, E.F., Protecting Turbomachinery from Self -Excited Rotor Whirl." 
Journal of Engineering for Power , Oct. 1965, pp338-344. 

2. Horlock, J.H. and Greitzer, E.M. , "Non-uniform Flows in Axial Compressors 
Due to Tip Clearance Variation," Proc. Instn. Mechanical Engrs., 

Vol. 197c , ppl73 

3. Urlichs, K. , Clearance Flow Generated Transverse Force at the Rotors of 
Thermo Turbomachines," Dissertation, Technical University of Munich, 
(1975). English translation in NASA TM- 7 7 292, Oct. 1983. 

4. Wohlrab , R., "Experimental Determination of Gap Flow-Conditioned Force at 
Turbine Stages and Their Effect on the Running Stability of Simple Rotors," 
Dissertation, Munich Institute of Technology (1975). English translation 
in NASA TM-77 293, Oct. 1983. 

5. Martinez-Sanchez, M, Greitzer, E. and Qiu, Y. , "Turbine Blade-Tip Clearance 
Excitation Forces," Report on NASA Contract Number NAS8-35018, June, 1985. 

6. Vance and Laudadio , "Experimental Measurement of Alford 1 s Force in Axial 
Flow Turbomachinery," Journal of Engineering for Gas Turbines and Power, 
1984, pp585 . 

7. Dunham and Came, "Improvements to the Ainley-Mathieson Method of Turbine 
Performance Prediction," ASME Transactions , July 1970. 

8. Lee, D.W.K., Martinez-Sanchez, M. and Czajkowski, E. , "The Prediction of 
Force Coefficients for Labyrinth Seals," 3rd Rotordynamic Instability 
Workshop, Texas A&M University, May 1984. 


List of Symbols 


c 



*xy 


P 


absolute velocity 
specific heat 
mean rotor diameter 
eccentricity 

blade force per unit circumferential length 
tangential force 
blade span 

excitation coefficient 
pressure 
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R blade radius 

t time; tip clearance 

T temperature 

T m mean torque on rotor 

w relative velocity 

u blade speed 

U s ideal blade force=( ideal turbine work)/(blade speed) 

a absolute flow angle 

$1 2 3 relative flow angle 

p ' ' Alford factor 

£ perturbation vorticity 

0 circumferential angle 

¥ stream function 

Q whirl speed 

p denity 

H efficiency 

Superscript 


perturbation 
mean value 


subscript 

y tangential direction 

x axial direction 

1 upstream of stator 

2 downstream of stator 

3 downstream of rotor 

0 stagnation 


rodial clearance// jlode height 



2 -6 % 


■HZ9I 

axial qap 

7 . 1 % 

.705 

.594 

.521 

*'% 

.716 

.570 

.440 


Table. I Efficiency for tip clearances 
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Figure 2. - Rotor under whirl. 


295 




tip clearance variation 
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Figure 4. - Axial velocity and local blade force distribution. 
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Figure 5. - Variation of excitation coefficient with tip clearance. 
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Figure 7. - Pressure perturbation distribution. 
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Figure 8. - Variation of excitation coefficient with tip clearance. 
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